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Abstract 

We discuss the possibilities of the simultaneous determination of the neutrino 
masses and the evolution of dark energy from future cosmological observations such 
as cosmic microwave background (CMB), large scale structure (LSS) and the cross 
correlation between them. Recently it has been discussed that there is a degeneracy 
between the neutrino masses and the equation of state for dark energy. It is also 
known that there are some degeneracies among the parameters describing the dark 
energy evolution. We discuss the implications of these on the cross correlation of 
CMB with LSS in some details. Then we consider to what extent we can determine 
the neutrino masses and the dark energy evolution using the expected data from 
CMB, LSS and their cross correlation. 
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1 Introduction 



Probing the masses of neutrinos has become one of the important targets in cosmology. 
Although neutrino oscillation experiments have well measured the mass differences as 
Am 2 2 ~7x 1CT 5 eV 2 from solar neutrino experiments [1-6] and Am| 3 ~ 2.6 x 1CT 3 eV 2 
from atmospheric neutrino experiments [7, 8] , they are insensitive to the absolute values 
of the neutrino masses. Experiments using kinematical probe such as tritium decay mea- 
surements can give an upper bound on an absolute neutrino mass, however cosmology can 
give a more stringent bound. From the analyses of recent cosmological observations such 
as cosmic microwave background (CMB), large scale structure (LSS) and so on, we can 
conservatively say the current bound on the sum of neutrino masses m v is around 2 eV 
(95% C.L.) [9-15]. 

Another important issue in cosmology is to understand the nature of dark energy. Al- 
most all current cosmological observations suggest that the present universe is dominated 
by an enigmatic component called dark energy. Although many models for dark energy 
have been proposed so far, we have not pinned down the model yet. However, by pa- 
rameterizing the dark energy with its equation of state wx, cosmological observations can 
give the constraints on wx- Assuming wx being constant, the current observations give 
Wx ^ _i [9-11,13,14,16-18]. Furthermore, even if we allow time-varying Wx and/or 
non-flat universe, it is shown that wx is constrained to be around a cosmological constant 
for some types of time variation [19]. 

It has been shown that the constraint on the neutrino masses from CMB, LSS, type 
la Supernovae (SNela) and so on can be weakened if the equation of state for dark energy 
is allowed to take a value other than —1 and vice versa [20], which means that there 
is a degeneracy between the neutrino masses and the equation of state for dark energy. 
In Ref. [20], the equation of state for dark energy is assumed to be constant in time, 
however, most models of dark energy proposed so far have a time-varying equation of 
state. Importantly it should be noticed that there are degeneracies among parameters 
which describe the time dependence of dark energy equation of state. For example, it was 
discussed that, for the models with a constant equation of state wx and a time- varying 
equation of state parametrized as wx = Wo + (1 — a)w±, there is a degeneracy between 
wq and w\ in the CMB power spectrum and matter power spectrum [21]. Since now we 
know that neutrino has a mass and the equation of state for dark energy can have the time 
dependence, we should take both of them into account when we consider the constraints 
from observations. However, considering the degeneracies discussed above, we can expect 
that we encounter unfortunate situations when we want to determine the neutrino masses 
and the evolution of dark energy simultaneously. 

In this paper, we discuss to what extent we can determine the neutrino masses and 
the evolution of dark energy, i.e., the time dependent equation of state, simultaneously 
from future cosmological observations. For this purpose, we make use of the Fisher matrix 
analysis using the expected data from future CMB and LSS observations. In addition 
to them, we also consider the cross correlation between CMB and LSS, which has been 
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attracting attention of the community in recent years. In particular, we investigate the 
effect of the neutrino masses on the cross correlation, which has not been attempted 
previously in the literature. As is well known, after the universe has been dominated by 
dark energy, the gravitational potential decays, which drives the so-called late integrated 
Sachs- Wolfe (ISW) effect. Since the time-evolving gravitational potential that drives the 
ISW effect may also affect large scale structure formation, the temperature fluctuation 
from the ISW effect and the distribution of galaxy are considered to be correlated [22-24] . 
Since there is no ISW effect in the universe dominated by matter, the detections of the 
cross correlation of the ISW effect with galaxy survey is a piece of physical evidence for 
the existence of dark energy. The detection of the cross correlation between the WMAP 
temperature fluctuation and several galaxy surveys have been reported in Refs. [25-35]. 
Many authors have discussed the cross correlation of CMB with LSS to investigate the 
properties of dark energy and some other issues [21,36-42]. In particular, it was shown 
that the cross correlation of CMB with galaxy survey can be a good probe for the evolution 
of dark energy [21]. 

The organization of this paper is as follows. In the next section, first we briefly discuss 
the degeneracy between the neutrino masses and the equation of state for dark energy 
and also that in models with the time-varying equation of state. Then, in section [3J we 
discuss the cross correlation of CMB with LSS and its spectra in scenarios with massive 
neutrinos and dark energy with some equations of state. In section HJ we study effects 
of massive neutrinos and dark energy on the suppression of growth of perturbation and 
the ISW effect and discuss how the masses of neutrinos and the equation of state for 
dark energy can affect the cross correlation spectrum. In section [5J we discuss a future 
constraint on the neutrino masses and the evolution of dark energy from CMB, LSS and 
the cross correlation of CMB and LSS. We give the summary of this paper in the final 
section. 

2 Degeneracies in neutrino and dark energy sectors 

In this section, we briefly discuss the degeneracy between the neutrino masses and the 
equation of state for dark energy wx, which was recently pointed out [20] and also the 
degeneracy in the parameters which describe the time dependent wx- 

First we discuss the degeneracy between the neutrino masses and wx- In Figs. [Hand 
El we show the CMB power spectra and matter power spectra respectively, for the cases 
with ACDM model and a model with massive neutrinos and the dark energy equation of 
state wx which is not equal to —1. For reference, we also plot the data from WMAP3 
[13] in Fig. Q] and those from SDSS [14] in Fig. H For the ACDM model, we take the 
cosmological parameters as Vt m h 2 = 0.13, Q^h 2 = 0.023, h = 0.735, r = 0.09 and n s = 
0.95 which are the mean value of the power-law ACDM models from WMAP3. Here 
Qi is the energy density for a component % {% can be b, m and u, which respectively 
stand for baryon, matter and neutrino) normalized by the critical energy density, h is the 
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Hubble parameter, r is the reionization optical depth and n s is the scalar spectral index 
of primordial fluctuation. In this paper, we assume that the tensor mode is negligible. 
For the matter power spectra, we have corrected for a scale-dependent biasing following 
the treatment of Ref. [14]. This correction affects only the scales with 0.07 /iMpc -1 . 
For the model with massive neutrinos and dark energy with wx 7^ — 1, we take the 
cosmo logical parameters as Q u h 2 = 0.01 and wx = —1.5. Here we assume degenerate 
masses for 3 neutrino flavors and a constant wx- The energy density and the masses of 
neutrinos are related by Q u h 2 = ^m„/93.2eV. The other cosmological parameters are 
taken as Vt m h 2 = 0.13, Q b h 2 = 0.023, h = 0.74, r = 0.087 and n s = 0.94. In general, 
dark energy component is characterized with its equation of state wx and speed of sound 
c 2 = Spx/Spx- Although the background evolution depends on wx only, the fluctuation 
of dark energy depends on both wx and c 2 . In this paper, we assume that c 2 = 1, which 
corresponds to the case with a scalar field dark energy, where (? s is defined at the rest frame 
of the dark energy. As seen from the figure, these two parameter sets give almost the same 
spectra, which means that the values of x 2 from CMB and LSS are almost the same. As 
shown in Ref. [20] , there is a degeneracy between m„ and wx in the direction to the region 
where wx < — 1. This degeneracy can be explained as follows. To fit a model with massive 
neutrinos to large scale structure data, large Q u h 2 can be compensated by increasing the 
energy density of matter. This cancellation seems not to be allowed because, assuming a 
flat universe and the cosmological constant for dark energy, SNela data constrain Vt m to be 
~ 0.3. However allowing the equation of state wx to vary, large values of Q m are consistent 
with observations of SNela with smaller wx [16]. This gives rise to the degeneracy in the 
matter power spectrum. As for CMB data, it is known that massive neutrinos affect the 
structure of acoustic peaks such as the position and the shape [12]. Increasing Q u h 2 , the 
peak position shifts to lower multipoles. However, by decreasing wx, the position of the 
peak shifts to higher multipoles, which cancels the shift caused by increasing the value of 
Q u h 2 . Although the shape of the peaks is also slightly modified by increasing fl u h 2 , by 
changing other cosmological parameters, we can fit such a model to CMB data well. Thus 
by varying wx, larger values of Q u h 2 become to be allowed by the data. 

Next we discuss the degeneracy in models with dark energy with a time- varying equa- 
tion of state. As mentioned above, most models of dark energy proposed so far such 
as quintessence have a time- dependent equation of state. Thus when we consider the 
constraints on the equation of state for dark energy wx phenomenologically, we should 
accommodate the time evolution of wx in some way. Although several ways have been 
proposed to include the time dependence of wx, here we assume a simple form as [11,43,44] 

wx(a) = w + (1 - a)wi. (2.1) 

In Figs. [T] and El the CMB and matter power spectra for the cases with wq = — 1 
and W\ = —1.2 are also plotted. Other cosmological parameters are taken to be the 
same as the case of a constant equation of state with massive neutrinos just above. As 
seen from the figure, the degeneracy in the CMB and LSS exists between the models 
with a constant wx and a time- varying wx- The degeneracy comes from the fact that 
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Figure 1: The CMB TT power spectrum for the ACDM model (red solid line), the model 
with wx = —1.5 and Q u h 2 = 0.01 (green dashed line) and the model with the time-varying 
equation of state w = —1 — 1.2(1 — a) and fl u h 2 = 0.01 (blue dotted line). Notice that 
three lines are almost indistinguishable. The data from WMAP [13] are also plotted. For 
the cosmological parameters taken here, see the text. 



both models give almost the same angular diameter distance from last scattering surface 
which determines the positions of the acoustic peaks in the CMB angular power spectrum. 
Furthermore, the fluctuation of dark energy is almost irrelevant to the structure of the 



acoustic peaksEj, hence they give almost an identical power spectrum provided the angular 
diameter distances from last scattering surface are the same. A possible difference between 
the models with constant and time dependent wx is the ISW effect which affects the low 
multipole region. In fact, there is a small difference in the region, however, it is not enough 
to differentiate them due to the cosmic variance in this case. 

Since there is a degeneracy between the neutrino masses and the constant equation of 
state for dark energy, we can expect that the situation becomes much worse when we also 
consider the time dependent dark energy equation of state. However, the possibility of 
differentiating the models with constant and time dependent wx using the cross correlation 
of CMB with LSS has been discussed [21]. Thus we may be able to break these degeneracies 
using the cross correlation. Before we discuss the future constraints on the masses of 
neutrino and dark energy evolution from CMB, LSS and their cross correlation, we study 

# 2 For the parameters assumed here, dark energy is negligible at the epoch of recombination. 
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Figure 2: Matter power spectrum for the same models given in Fig. [TJ The data from 
SDSS [14] are also plotted. 



the effects of them on the cross correlation of CMB with LSS in the next section. 



3 Cross correlation of CMB with galaxy 

In this section, first we briefly review the formulation of the cross correlation of CMB with 
galaxies following Ref. [36]. Detailed descriptions of this issue can be found in Refs. [23, 
24,36]. 

The two point correlation function between the ISW signal and galaxy overdensity is 
defined as 

a Isw ^(e) = (A IS w(ni)5 g (n 2 )}, (3.2) 

where the angular brackets represent the ensemble average and cos# = hi ■ h 2 . A ISW is 
the temperature fluctuation from the ISW effect in the direction hi which is written as 

Aisw(ni) = = -2 J e-^&(hi, r))d V , (3.3) 

where rj is the conformal time and e _T ™ is the visibility function. Here a prime denotes 
the derivative with respect to the conformal time. $ is the gravitational potential which 
appears in the metric perturbation in the conformal Newtonian gauge as 

ds 2 = -a 2 (l + 2§)drj 2 + a 2 (l + 2^)dx 2 . (3.4) 
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Unless otherwise stated, we work in the conformal Newtonian gauge in this paper. 
5 g (n 2 ) is the overdensity of galaxies in the direction n 2 . 

W = ^4^, (3.5) 
n g 

where n g is the mean number density of galaxies. It is assumed that the galaxy number 
overdensity traces the density fluctuation of matter as 

5 g = bS m , (3.6) 

where b represents the galaxy bias. Thus 5 g can be written as 

5 g (h 2 ) = b J (f)(z)z'S m (n 2 ,v)dv, (3-7) 

where <f>(z) is the selection function of a given galaxy survey. 

As usual, it is convenient to decompose C ISW ~ g (9) with the Legendre polynomials 
Pi(cos6) as 

00 9/ 4- 1 

c isw- g{e) = J- ___tic; SW - g P,(cos0). (3. 8 ) 

1=2 

We subtracted the monopole and dipole contributions by construction. With the above 
definition, the cross-correlation power spectrum is given by 

c™-*=4*±J^irmw, 0.9) 

where is the primordial power spectrum. Ij sw (k) and If(k) are defined as 

lF W (k) = ~2 J e-^H'Mmm ~ V)]dv, (3.10) 

If(k) = b J (f)(z)z'5^ n (r])ji[k(r] - rj^drj, (3.11) 

where the ji is the spherical Bessel function and rj is the conformal time at the present 
epoch. $ fc and are the Fourier component of the gravitational potential and density 
fluctuation of matter which can be written as 

/H 3 k 
j^y 3 $k(ny k - h \ (3.12) 

/d 3 k 
J^yA(n)e^, (3.13) 

where r = r/Q — r]. For the later use, we also write down here the galaxy-galaxy correlation 
function which is given as 

w (U) =47r ^ J d ^ n lf {k )lf{kl (3.14) 




where indices represent the redshift bins for the galaxy selection functions. 

Now we show the cross correlation spectrum in models introduced in the previous 
section. We computed the cross correlation spectra using the modified version of CMB- 
FAST [45]. In Fig. [3l we plot the cross correlation power spectrum for the same models 
as those in Fig. [TJ Here we use the normalized Gaussian selection function with the peak 
at redshift z = 0.5 and the variance a z = 0.07. The other cosmo logical parameters are 
the same as those in Fig. [TJ In the figure, the scalar amplitudes are normalized to the 
WMAP best fit values and the bias factors are determined by the fit to SDSS data. We 
will investigate in Sec. [5] whether or not the difference seen in Fig. [3] can be probed by 
future cross-correlation data. Before that, in the next section, we investigate the cross 
correlation in detail, discussing the suppression of matter fluctuation and the ISW effect 
in models with massive neutrinos and dark energy with some equations of state. 

4 Suppression of perturbation growth and the ISW 
effect 

In this section, we discuss the effects of massive neutrinos and dark energy with some 
equations of state on the suppression of fluctuation growth and the ISW effect. Below, the 
evolution of perturbation variables are calculated at the mode k = O.Olh Mpc -1 . Then we 
study the cross correlation spectrum in the models. 
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Figure 4: Suppression of fluctuation growth in models with dark energy equations of state 
w x = —1 (red solid line), wx = —0.8 (green short-dashed line), wx = —1.5 (blue dotted 
line) and w x = —1 — 1.2(1 — a) (purple long-dashed line). In these cases, we assume 
that VLyh 2 = 0. The case with massive neutrinos Q u h 2 = 0.01 and wx = —1 (light blue 
dash-dotted line) is also plotted. 



First we consider the effects of massive neutrinos on the suppression of fluctuation 
growth. Since massive neutrinos represent a smooth gravitationally stable component on 
small scales, they suppress the growth of matter fluctuation [46]. In linear regime, a small 
initial fluctuation <5 grows as 5(a) = D(a)5o where D(a) is the growth factor. During 
matter-dominated universe, D(a) grows as D(a) oc a. An analytic formula for the growth 
factor is given in Ref. [46]. In Fig. HI we show S m /a for the model with Q u h 2 = 0.01 and 
Wx = —1 as a function of redshift. Other cosmological parameters are taken to be the 
same as those of the ACDM model in Fig. [TJ We can see the suppression of growth of 
matter fluctuation in models with massive neutrinos. The growth of matter fluctuation 
is more suppressed by increasing the masses of neutrinos. To see the ISW effect in the 
model, we also plot d^/drj as a function of redshift in Fig. [5j As is clear from the figure, 
massive neutrinos have little effect on the evolution of the gravitational potential. Thus 
the mass of neutrinos is almost irrelevant to the late ISW effect, which means that we can 
expect that the neutrino masses have almost no effect on the cross correlation of CMB 
with LSS except for the overall suppression of matter fluctuation. Notice that, since we 
do not have much knowledge of the bias factor, we might not see such suppression. 

Next we discuss the effects of dark energy on the suppression of fluctuation growth 
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Figure 5: The derivative of the gravitational potential with respect to the conformal time 
for the same models given in Fig. HI 



and the ISW effect. First we will discuss the suppression of fluctuation growth. The 
growth of matter fluctuation can be described by the following equation [47], including 
the fluctuation of dark energy 

Sc " + m' c - ^fn c 6 c = ^fn x s x , (4.15) 

where TC = a' / a is the conformal Hubble parameter, 5x is the density perturbation of dark 
energy and fli(z) = pi (z) / ptot&\(z) is the energy density of a component i normalized by 
the total energy density at redshift z. 

In Fig. HI the fluctuations of matter for the cases with wx = —1.5,-1,-0.8 and 
wx = — 1 — 1.2(1 — a) are shown. The matter fluctuation is more suppressed as increasing 
the value of wx- This is due to the change of the background evolution by the domination 
of dark energy. The fluctuation of dark energy also affects the fluctuation of matter via 
gravity. 

Another important effect of dark energy is the late ISW effect. As is welhknown, when 
the universe is dominated by dark energy, the gravitational potential decays, which drives 
the late time ISW effect. The ISW effect can be written as 

$ ' ~ ^ [H 2 (tt m 5 m + Q x Sx)} ■ (4.16) 
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In Fig. [5j the evolutions of $' as a function of redshift for several values of wx are shown. 
Increasing the values of wx, the ISW effect is more enhanced at high redshift region, 
which leads to enhanced cross correlation power spectrsF 3 !. Notice that, at recent epoch, 
the values of $' decreases as the values of wx increases. Thus the effects of wx can depend 
on the values of wx itself and also on the redshift. 

Now we are in a position to discuss the cross correlation power spectrum for models 
with massive neutrinos and dark energy with some equations of state. In Fig. [BJ we show 
the cross correlation power spectra for the cases with the constant equations of state 
(wx = —0.8, —1 and —1.5) and the time- varying equation of state wx = — 1 — 1.2(1 — a). 
The case with Q u h 2 = 0.01 and wx = —1 is also plotted. We use the same galaxy selection 
function as that used in Fig. [31 Other cosmological parameters are taken to be the same 
as those of the ACDM model in Fig. [IJ Since the masses of neutrinos suppress the growth 
of fluctuation, we can see the power spectrum of the cross correlation is also suppressed. 
Although the growth of matter fluctuation can also be suppressed in the dark energy 
dominated universe, effects of dark energy also come from the ISW effect. As the equation 
of state of the dark energy increases, the epoch where the energy density of dark energy 
dominates the universe becomes earlier than the case with the cosmological constant, 
which affects the evolution of the gravitational potential more at high redshift. Thus it can 
produce the large ISW signals. However, as mentioned above, the size of the ISW effect also 
depends on the epoch we consider. Thus the size of the cross correlation power spectrum 
depends on the values of wx, its time dependence (i.e., w\ in our parameterization) but 
also on the selection function. 



5 Future constraints 

Now we discuss attainable constraints on neutrino masses and the evolution of dark energy 
from future experiments of CMB, LSS and the cross correlation between them. We adopt 
the usual Fisher matrix method [48] for this purpose. The Fisher matrix is defined as 

where C is the probability of observing a data set {xi,X2, ■■■} for a given cosmological 
parameter set {Ai, A2, ••••}• Since there have been many works on the Fisher matrix analysis 
using future CMB and LSS data, we refer Refs. [49-51] for the details. 

As for the cross correlation, here we briefly summarize the formulation following 
Ref. [21]. To forecast the errors in cosmological parameters, we use the expected CMB 

# 3 In fact, if we consider a smooth dark energy, the tendency is different. When we assume a smooth 
dark energy, the effect of dark energy on the ISW effect comes from the modification of the background 
evolution alone. Thus in this case, increasing wx makes the epoch when the universe begins to accelerate 
later, then the ISW effect becomes less. 
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Figure 6: Cross-correlation power spectra for the same models given in Fig. HI 



data from PLANCK [52] and the expected galaxy survey by Large Synoptic Survey Tele- 
scope (LSST) [53]. For PLANCK, we use the experimental specification which can be 
found in Ref. [54]. For the distribution of galaxies from a given survey, we assume that 
the total galaxy number is given by [21,38], 



(5.18) 



where z n is a parameters which gives a median redshift of the survey. The galaxies can be 
subdivided into multiple bins as 



< tal w = E^)> 



(5.19) 



through photometric redshift. To approximate the redshift binning, we assume that 
photometric redshift estimates are Gaussian distributed with the variance of a(z) = 
o"max(l + z)/(l + ^max) [21]. Then the photometric redshift distributions are given by 



erfc 



Zi-l ~ z 



y/2, 



viz 



erfc 



Zi - z 



y/2, 



(TiZ 



(5.20) 



where erfc is the complementary error function. For the LSST survey, we assume the 
galaxy survey sky coverage fraction as / s k y = 0.5 and 10 photometric redshift bins out to 
z ~ 3. The total galaxy number density is assumed to be 70 gal/arcmin 2 . 
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For the CMB-galaxy cross correlation, the Fisher matrix can be written as [21] 



FT'' = f*, E< 2 ' + d £ -§—(c ^i-^-. (5.2i) 

where the covariance matrix (Covj)isw-g is given by 

l^ OV UlSW-g - °Z U l + °2(i) • O^J 

Here is the observed spectrum which includes the noise 

-l 



cr = cr + £ 



e 



ff C g FWHM,C \ ,/i/-L)«Mvn:.j.r/Sl.i2 
?CMB 



(5.23) 



where ^ c represents the summation over channels, oq and #fwhm,c are the sensitivity 
per pixel and the width of the beam of a channel for a given measurement, respectively. 
Tcmb is the temperature of CMB. cDj is the observed auto-correlation which is the sum 
of the signal and the Poisson noise 

Co^ = uf' 3) + Nj l ' J \ (5.24) 
where N^ 1 '^ is the Poisson noise which is uncorrelated between bins, thus 



N^ j) = (5.25) 



n\ is the galaxy number per solid angle in the i-th bin. For the future constraints from 
CMB and LSS, we followed the method presented in Ref. [54]. 

Now we discuss attainable constraints from future experiments. As the fiducial model, 
we take the cosmological parameters as Qj,h 2 = 0.024, Q m h 2 = 0.14, Qx — 0.73, r = 0.166 
and n s = 0.99. For the neutrino masses, we take m u = 0.3 eV with the three degenerate 
masses. We checked that even if we assume the normal hierarchy with smaller neutrino 
masses, our conclusion does not change much. For the dark energy equation of state, we 
take Wx = —1 for the constant case and wx = — 1 + 0.3(1 — a) for the time dependent 
case. In Fig. [3, we plot the expected constraints on the neutrino masses and dark energy 
equation of state for the cases with the constant and the time dependent wx on the left 
and right panels respectively. When we obtain the constraints on the Q u h 2 vs. wq plane, 
we marginalized over other cosmological parameters. Especially, for the latter case, we also 
marginalized over W\ to obtain the constraint. In the figure, la CL. contours are shown 
for the future data from CMB only, CMB+cross correlation, CMB+LSS, and all combined. 
As we can see, for the case with constant %, the cross correlation observations do not 
help much to constrain the neutrino masses and dark energy equation of state. However, 
for the case with the time-evolving dark energy equation of state, the cross correlation can 
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Figure 7: Expected la contours on the VL v h 2 vs. Wq plane. The cases with the constant 
Wx (left) and time-varying wx (right) are shown. We marginalize W\ for the latter case. 
Contours which correspond to the attainable future constraints from CMB only (red solid 
line), CMB+cross correlation (blue dotted line), CMB+LSS (green dashed line) and all 
combined (purple dash-dotted line) are shown. 



be useful to determine the time evolution of wx- As for the neutrino masses, the constraint 
on m u becomes weaker just a little in the case with the time evolving Wx- In both cases, 
the inclusion of LSS data can help to constrain the neutrino masses. As mentioned above, 
even if we assume the normal hierarchy for the neutrino masses, the shape and the size 
of the contours scarcely change. In fact, even if we marginalize over u>i, the constraint on 
the neutrino masses is not so weakened. This means that the time dependence of wx (i.e., 
Wi in our parameterization) and m v are not strongly correlated. 

Next we discuss the expected constraints on the dark energy evolutions, i.e., w and 
W\. In Fig. [HI we plot the future constraints on the time evolution of dark energy, fixing 
the neutrino masses (left panel) and marginalizing over the neutrino masses (right panel). 
As we can see, when we marginalize the masses of neutrinos, the uncertainty of Wq and 
Wi becomes larger, in particular, in the direction parallel to the Wq axis. In fact, the 
constraint on W\ is not weakened much compared to that on Wq. As mentioned above, 
this is because the masses of neutrino and W\ are not so degenerate, which means that the 
time dependence of the equation of state for dark energy is not strongly correlated with 
the neutrino masses. It should also be noticed that the cross correlation of CMB with 
LSS can help to some extent to determine the evolution of dark energy even when we have 
uncertainties in the absolute values of the neutrino masses. 
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Figure 8: Expected la contours on the Wo vs. w\ plane. On the left panel, we fixed 
the neutrino masses to the fiducial value. On the right panel, we marginalized over the 
neutrino masses. Contours which correspond to the attainable future constraints from 
CMB only (red solid line), CMB+cross correlation (blue dotted line), CMB+LSS (green 
dashed line) and all combined (purple dash-dotted line) are shown. 



6 Summary 

We discussed the simultaneous determination of the neutrino masses and the evolution 
of dark energy from future cosmological observations. After presenting the degeneracies 
which exist between the neutrino masses and dark energy sectors, we considered the pos- 
sibilities of breaking the degeneracies using the cross correlation of CMB and LSS. We 
studied the cross correlation power spectrum in those models in some details. We also 
discussed the suppression of perturbation growth and the ISW effect in the models. Then 
we presented the future constraints on the neutrino masses and the dark energy equation 
of state from the Fisher matrix analysis using the expected data from CMB, LSS and their 
cross correlation. First we showed the results on the Q u h 2 vs. Wq plane for the cases with a 
constant equation of state and a time-varying equation of state. We can see that the cross 
correlation of CMB with LSS can help to determine the time dependence of the equation 
of state to some extent. The constraint on the neutrino masses is not affected much since 
there is not strong correlation between the neutrino masses and w± which represents the 
time dependence of dark energy equation of state in our parameterization. Then we also 
discussed the constraints on the time evolution of equation of state, fixing and marginal- 
izing the neutrino masses. When we marginalize the neutrino masses, the constraints on 
wq and Wi are weakened. However, we showed that the cross correlation observation can 
be useful to determine the time evolution of the equation of state even if we do not have 
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much information on the masses of neutrinos. 

Note added: Recently the paper [55] appeared on the arXiv, which discusses the effect 
of the neutrino masses on the cross correlation. Also, Ref. [56] appeared more recently 
which has some overlaps with our analysis. Ref. [56] gives a qualitatively similar result to 
ours that the effects of the neutrino masses on the cross correlation are relatively small. 

Acknowledgment: We acknowledge the use of CMBFAST [45] package for our numerical 
calculations. T.T. would like to thank the Japan Society for Promotion of Science for 
financial support. 
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